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Abstract. We give asymptotic formulae for random matrix averages of derivatives of char- 
acteristic polynomials over the groups USp(2iV), SO(27V) and 0~(2N). These averages are 
used to predict the asymptotic formulae for moments of derivatives of L-functions which 
arise in number theory. Each formula gives the leading constant of the asymptotic in terms 
of determinants of hypergeometric functions. We find a differential recurrence relation be- 
tween these determinants which allows the rapid computation of the [k + l)-st constant in 
terms of the fc-th and (k — l)-st. This recurrence is reminiscent of a Toda lattice equation 
arising in the theory of T-functions associated with Painleve differential equations. 



For over 50 years, mathematicians and physicists have used random matrix theory to 
study a wide-ranging, growing list of probabilistic phenomena. Particularly surprising are 
its applications in number theory, where random matrices model the distribution of nontrivial 
zeros of the Riemann zeta function. Random matrix theory now provides far-reaching and 
widely believed conjectures for many questions in the analytic theory of L-functions. 

Katz and Sarnak [KaSaJ give evidence that every family of L-functions falls into one of four 
symmetry types: unitary U(iV), unitary symplectic USp(2iV), even orthogonal SO(2A r ) and 
odd orthogonal 0~(2N). These symmetry types govern the distribution of zeroes and special 
values in families. Using random matrix models, Keating and Snaith |KeSn] and Conrey, 
Farmer, Keating, Rubinstein and Snaith |CFKRS0~3] have produced deep conjectures for 
estimating the integral moments of central values in families of L-functions. 

The derivatives of L-functions are also of great interest, and are the subject of this paper. 
A motivational example is Speiser's theorem, which asserts that the Riemann hypothesis is 
equivalent to the nonexistence of nonreal zeros of the derivative of the Riemann zeta function 
to the left of the critical line, see e.g. [SouJ. Moreover, the derivatives of L-functions control 
the order of vanishing at the central point, which encodes important arithmetic and geometric 
information. For example, according to the Birch and Swinnerton-Dyer Conjecture, the 
order of vanishing of the L-function of an elliptic curve over the rationals coincides with the 
arithmetic rank of the curve. 

An L-function is modeled by the characteristic polynomial of a random matrix A. Here 
we compute 
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where G denotes USp, SO, or O , and dA is the Haar measure on G. As N — > oo, this 
models the kth moment of L^ m ^(l/2) in a family of symmetry type G. 



One can find the moments of A^(l) by differentiating the corresponding shifted moment 
formulae, which are computed in |CFKRS03] . Conrey, Rubinstein and Snaith in |CRS] 
develop a faster method to compute the relevant averages in the unitary case: 



/ \A' A (l)\ 2k dA = b k N k2+2k + (N k 

JV(N) ^ 



T k 2 +2k-l 

'V(N) 

This leading constant bk is the same "geometric constant" appearing in conjectures for the 
asymptotic estimate of the fc-th moment of ('{s) in t aspect. Conrey, Rubinstein and Snaith 
describe bk in terms of a k x k determinant of /-Bessel functions, and are able to compute 
bk numerically for k < 15. 

Forrester and Witte in |FW06tlFW02] find a surprising expression for these determinants 
of /-Bessel functions in terms of solutions to Painleve III' differential equations. Explicitly, 
the formula in |CRSj is as follows: 

k /,\ / j x h+k 



t=0 



where Ie(x) is the modified Bessel function of the first kind. One then defines 

T k {t) := 2- fc ( fc - 1 >r fc2 / 2 det (h^-^Vi)) . 

k x k V J 

Forrester and Witte |FW06llFW02] find that this r k (t) (denoted r[k](t) in jFW02l section 
4]) is in fact the Okamoto r function associated with the Painleve III' differential equation: 

(ty"f + y'(Ay'-l)(y-ty')-^k 2 = 0. 

This nonlinear second order differential equation has a solution with certain boundary data 
(sec [FW02J) given in terms of Tjt(t) by the formula 



y = amAt) = -tj f \og ( e "* /4 ^(i)) 



Specifying boundary conditions, one can quickly compute (Tm,k(t) from the differential equa- 
tion and recover T&(t) via the equation 



r k (t) = exp (- J 



(\ 7 2 \ 

s)+k - -) — 
4 s 



This expression allows a much faster computation of the constants bk- 

The goal of the present paper is to extend these results to the other symmetry types 

relevant to L-functions: USp(2iV), SO(2i V) and 0~(2N). 

Employing similar techniques to [CRSj . we obtain analogous results, where the role of the 

/-Bessel functions above is here played by hypergeometric functions, 

1 f e*"^^ 
9m(u) = — <b —^rrdw 



2ni J lw]=1 w m+1 



r(m + iy * V 2 
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for 11 6 C and m G Z. For negative m, interpret the above expression as the limit. The role 
of the r-function is played by 

T k /{u) := det (g 2 i- j+i {u)) , (1) 

kxk 

for k > 0, £ G Z and u G C, where, here and in the following, the indices i and j of the matrix 
in the determinant range from 1 to k. In the context of Theorems [IJ [2] and [3] below, the i 
appearing here takes values 0, —1, respectively. We now state theorems for the symplectic, 
special orthogonal, and negative orthogonal cases: 

Theorem 1. We have 

M fc (USp(2iV),2) = 6 fc (USp(2iV),2) ■ (2N)^ + O(N^) (2) 



where 



6 fc (USp(2iV),2) = 2-"^^ (e"T M (2u))| u= 



Theorem 2. We have 

M fc (SO(2iV),2) = 6 fc (SO(2iV),2) • (2N)^~ + 0(N^~) (3) 



where 



6 fc (SO(2iV),2) = 2-^^ (e«r fe> _ 1 (2 U ))| u= 



Theorem 3. We have 

M fc ((T(2A0,3) = 6 fc (0"(2iV),3) • (2N)^~ + 0(N~^) (4) 

where 

6*(CT(2A0,3) = 3 ■ 2^^ (e«r fc ,o(2n))| u=0 . 

Note that we find above that b k (0~(2N), 2) = 3 • 2 k ■ & fc (USp(2JV), 2). 

We are naturally led to consider the second and third derivatives of characteristic polyno- 
mials in the above theorems instead of the first derivative due to root number considerations. 
Indeed, if A is a unitary matrix, the characteristic polynomial satisfies the functional equa- 
tion 

A A (s) = (sfidetAy'AAis- 1 ), 



where f(s) = f(s). When A is in USp(2iV) or SO(2iV) then det A is constantly equal to +1 
and one has a simple expression for A^(l) in terms of Aa(1), which can be used to compute 
the moments of the derivative via partial integration. Thus the moments of A'^(l) give the 
next novel information. When A G _ (2AT) we have det A = —1, and thus A^(l) = 0. 
In this case A^(l) plays the role which A^l) plays in the other families, A^(l) has a 
simple expression in terms of A^(l), and therefore one considers moments of A^'(l). The 
same reasoning carries over to the families of L-functions having each of the aforementioned 
symmetry types. 

The method of proof of Theorems [TJ [2] and [3] can be generalized to higher-order derivatives 
easily. It suffices to expand the binomial in Lemma HJ and otherwise proceed as in the given 
proofs of Theorems [TJ [2] and [3J 



We can use Theorems [TJ E] and [3] to give conjectures for moments of derivatives of L- 
functions at s = 1/2 with the above symmetry types. For example, the quadratic Dirichlet 
//-functions ordered by conductor form a symplectic family and thus we make the following 
conjecture: 

Conjecture 1. Let T>(X) := {\d\ < X, d fundamental discriminant}, and L(s,Xd) denote 
the quadratic Dirichlet L-function of fundamental discriminant d. The average value of the 
second derivative of quadratic Dirichlet L-functions at the central point is 



L"(l/2, Xd ) k ~ a k ■ 6 fe (USp(2iV),2) • (logX)^ 



\V(X)\ 

Here a k is a well understood arithmetic constant depending on the functional equations 
in the family: 

(l-l/p)W^(l-l/^ + (l + l/^)-* 1 

ak ~ 11 1 + 1^ I 2 + p 

p prime ' N 

(see, for example (CFKRS05J 1.3.5]). This is the same arithmetic constant appearing in 
moment conjectures for L(l/2,Xd) without the derivative. It is not predicted by random 
matrix calculations. 

Our Theorems [H, [2] and [3] along with the results of [CRS] allow similar conjectures to be 
made for any family of L-functions. 

Under closer examination, the determinants Tk,t{u) defined by ([1]) and appearing in The- 
orems [U [2] and [3] exhibit a surprisingly rich structure. Our Theorem H] is a differential 
recurrence relation which allows much faster computation of the constants &&(G(iV), m): 

Theorem 4. Let k e Z >0; t G Z. Then 

T k+ i,i(u)Tk-iAu) = 2 (uT k ,z(u)T k " e (u) + TkAu)T k \ e (u) - u (T k \ e (u)) 2 ^j . (5) 

This recurrence relation closely resembles a Toda lattice equation for the Okamoto r- 
function associated with a Painleve differential equation, see [Okaj Theorem 2]. Such Toda 
lattice equations are at the heart of the r-function theory of Painleve equations, and are used 
by Forrester and Witte [F W0 6, FW0 2) to connect determinants the of i-Bessel functions 
found by [CRSJ to the Painleve III' equation. It would be very interesting to determine 
whether or not there exists a differential equation arising from our formula (jSJ) which plays 
the role for symplectic and orthogonal types that Painleve III' plays for unitary symmetry. 

Ultimately, one hopes to obtain formulae for the complex moments of characteristic poly- 
nomials in the various symmetry types. In the case of the undifferentiated moment conjec- 
tures, it has been found that the geometric constants g k can be expressed in a simple form 
in terms of Barnes G-functions, which are well-defined for complex values of k, see [CF] . A 
project for the future would be to see if there exists a similar expression for the geometric 
coefficients bk studied in this paper for moments of derivatives of L-functions. 

By computing the expressions in Theorems [lj [2] and [3] directly, one can obtain the values of 
bk up to k ~ 10. By using Theorem H] we do much better: running SAGE for about an hour 
on a machine with 4 gigabytes of RAM we computed the first 200 values of bk(USp(2N), 2). 
In section H] we give a table with the first 10 values of bk for each symmetry type. To give 
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an example, we have 

MUSp(2A0) = 



47 ■ 1553 ■ 1787 ■ 73709 ■ 152825093 

2 62 . 334 . 517 . 710 . U 5 . 135 . 17 4 . 19 3 . 23 2 • 29 • 31 • 37' 



, , nn , 1 « 25171 • 7695491 • 57668937071891 

6io(SO(2iV)) 



6 10 (O"(2iV)) 



2 45 . 329 . 515 . 79 . n s . 13 5 . 17 3 . 19 3 . 232 . 29 • 31 • 37 
47 • 1553 • 1787 • 73709 • 152825093 



252 . 333 . 517 . 710 . n s . 135 . 17 4 . ig3 . 23 2 ■ 29 ■ 31 • 37' 
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2. Some lemmas 

We recall the definitions of the relevant spaces of matrices: USp(2iV) is the subgroup of 

2N x 2N unitary matrices M with M* ( ^ ^ ^ M = ^ ^ ) ' w ^ ere ^ denotes ^ ne 

N x N identity matrix. SO(2iV) and 0~(2iV) are the subsets (O - is not a subgroup) of 
orthogonal matrices with determinant 1 and —1, respectively. These compact spaces admit 
Haar measures which we normalize so that the volume of each space is 1. A matrix in 
USp(2iV) or SO(2iV) has characteristic polynomial of the form A(x) = Yl n =i(^ ~ e ie ™x)(l — 
e~ "s), and a matrix in 0~(2N) has characteristic polynomial of the form A(x) = (1 — 
x){l + x) Il^i 1 ! 1 - e Wn x)(l - e- w -x), with 9 n e E. 
The shifted moments of A are defined as follows: 



I(G{2N);z 1 ,...,z k ):= / A( Zl ) ■ ■ ■ A(z k ) dA. 

JG(2N) 

Conrey, Farmer, Keating, Rubinstein, and Snaith |CFKRS03] use Weyl integration formula 
(see, for example, Theorem 8.60 of [Kn]) to compute the following shifted moment formulae, 
which are the starting point for our work (note that we have corrected a typo in [CFKRS03~1 
4.9]): 

Lemma 1. [CFKRS03, 3.36] Assume that a 1; . . . , are complex numbers with |cK f | < 1 for 
i = 1, . . . k. Then 

J(USp(2iV);e- Ql ,...,e- afc ) = 

( _ 1)fc (*-i)/2 2 * r r YlxKi^f ~ UU W - 



{2m) k k\ J {ml=1 J ]Wk]=1 Ui<ij<k( w j ~ a i) 



3 

x 



x e iVE^i Yl (1 - e-^-^y 1 d Wl --- dw k . 



l<m<e<k 
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Lemma 2. [CFKRS03, 4-43] Assume that a±, . . . , are complex numbers with \aii\ < 1 for 
i = 1, . . . k. Then 

J(S0(2iV);e- ai ,...,e- Qfc ) = 

x 



{2ni) k k\ J {wil=1 J ]Wk]=1 Ui<i,j< k ( w ] ~ a i) 

x e ^E* =1 K+«,) ] [ (i _ e -«m-«*)-i dwi . . . rfWfe . 

l<m<^<fe 

Lemma 3. \CFKRS03i 4-9] Assume that a>i, . . . , a k are complex numbers with |c^| < 1 for 
i = 1, . . . k. Then 

I(0-(2N);e- a \...,e~ ak ) = 

( _ 1)fc(fc -i)/2 2fc , , n^^K 2 - ^ 2 ) 2 n -=i a 



3 

X 



\<m<i<k 

In fact, we will use approximate versions of the lemmas above, which follow immediately 
from the fact that (1 — e~ x )~ 1 = x^ 1 + 0(1). To simplify notation, we denote Vandermonde 
determinants as follows: 

A(w) := det(^ _1 ) = I I (Wi - wj) 

kxk - 1 - -•- 

and we write w 2 = {wf)i<i<k for any w = (Wi)i<i<k G C k . Then we have: 

Corollary 1. Assume that ati,...,ack are complex numbers such that \atj\ <C ■h for j = 
1 . . . k. Then 

J(USp(2iV);e- Ql ,...,e- Qfe ) = 

(_l)fc(fc-i)/2 ft r A{w)A{w 2 )e N ^^- a ^ 



> . w 7 t i t r i-r ■> ■■>> -dw 1 ...dw k \(1 + 0(N 
(2mfk\ \J j ll,. ; f<7 -of) ) K 

Corollary 2. Assume that ai,...,ak are complex numbers such that \aj\ <C for j 
1 . . . k. Then 

I(S0(2N);e- ai ,...,e- ak ) = 

= W n~R^?) d Wl ... dWk j(l + 0(N )). 

Corollary 3. Assume that ax,...,(Xk are complex numbers such that \aj\ <C for j 
1 . . . k. Then 

J(CT(2iV);e- Ql ,...,e- afc ) = 

{ _ l)k (k-i),2 2 k ( j. j- A(w)A(w 2 )(U j a j )e N ^^ 

' ' >?> _ dw\ . . . dw k (1 + U[N )). 



6 



The following formula will appear in the proofs of Theorems [T] and |2j below. We use it 
only when m — 2, but the general statement is a starting point for computing moments of 
the mth derivative. 



Lemma 4. For m > 0, we have 



d" 



d 1 - 



da™ da m ] [ 



l<i,j<k\ j 



w 2 — a?) 



ai=---=a k =0 



e=o 



E 7 (-")-' E n 



by 



h+-+i k =e, j=i w j 



ij even 



Proof. We have 
d m 



-Na 



dam Ui<3< k (tf-a 2 ) 



e=o 



e 7 ^ n 



L c?a l J (w| — a 2 ) 



and 



da 9 \w 2 -q; 2 



a=0 



if g is odd 

q\/w q+2 if g is even. 



Therefore, 



-Na 



da m H 



= E 

<*=0 <=0 



n? 



m—l„—Na 



_ N ^m-i e 



e n 



ij+2 



h+-+i k =£, j=l W j 



and the lemma follows. 



□ 



Our proofs will also rely upon Vandermonde determinants of differential operators: 



d 

dx 



n 

l<i<j<k 



d? 



-i 



d d . 
— = det I — _ . 



We give two lemmas on computing with these-Lemma [5] is a direct consequence of the 
definition, but we prove Lemma [6] in detail. 

Lemma 5. Let fi(x), . . . , fk(x) be k — 1 times differentiable. Then 



d 



1=1 



(xi) = det /, 

kxk 



(i-i), 



Lemma 6. Let f(x,y) be k — 1 tames differentiable in x and y. Then 

k 



A 



d 

dx 



A 



d 
dy 



\\f (■'■,!!, 



i=i 



xi=...=x k =X, 



fcxfc \ dX % ~ l dY3 



Proof. By Lemma [5] we have 

/ d \ k ( d?" 1 \ k d^^ 1 

a [fa) n H x *> v*) = ti \ahF~ if ^ yi) ) = s sign(/i) n ^FT /(Xi ' m) 

where the sum runs over the permutations fi G Sfc. Applying Lemma [5] again, we find 
d \ . / d \ A „, , , w / d\ A 



/ ^(j)+i-2 



and so 



M£) A G£)n/<*.K 



Now, we may rearrange the rows of the matrix 



to obtain 



doing so cancels out the sign(/x) attached to the determinant, and we reach the desired 
formula. □ 

Our final lemma is a recursion for determinants discovered by Lewis Carroll. It will be 
used in the proof of Theorem HI 

Lemma 7. Let A be an k x k matrix, and let ^4(^''"'^ r ) denote the matrix A with rows 
ai, . . . , a r and the columns b\, . . . , b s removed. Then 

(1(1 .4 ( I • ( Let A I { ] - det A r ) • det A ( { ] = det A ■ det A ' ' ' 7 



J/ \3J W \M 

3. Proof of the theorems 



A simple calculation shows that 

d m d m 



da? da™ v v ' n ai =...a k =o 



I(G(2N);e- a \...,e- ak 

n 

/* / " L f "\ \ 

mk 



:-iy 



772 I 

where < . > denotes a Stirling number of the second kind. It follows that 

U J 

M fc (G(2A0, m) = ^ • • • ^ J(G(2iV); e~-, . . . , e'^^^ (1 + O(N^)). 
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Thus we may proceed towards Theorems HJ El and [3] by differentiating the formulae found 
in Corollaries HJ El and |3l respectively The asymptotics remain valid after differentiating 
because they are uniform in aj. 



Proof of TheoremUl From Corollary [T] and the above argument, we know that the kth mo- 
ment of A" is asymptotically 

M fc (USp(2iV), 2) = ^i- M fc (USp(2iV), 2)(1 + O^- 1 )) 

/e! 

where 



M fc (USp(2iV),2) 



d 



2/,- 



A(w)A(w 2 )e Ar ^' ( ^-^ ) 



da\--- da\ (2mf J J [LjK - «?! 
We apply Lemma H] with m = 2, finding 



du>i . . . dtOfc 



cki •••=«& =Q 



d 



2k 



daf • • • da\ \\ x 



<i,j<k\ w j 



a?) 



a 1 =—=a k =0 



d£* 



k \ / fc^ 



1 W 7 



1 W ? 



t=0 



(6) 



This allows us to write 



jfe tTV 2 

M fc (USp(2iV),2) = — 



dt k (2iri) k 



j-.-j A{w)A{w 2 ) exp (^T Nwj + ^ j ^ 



du>fc 



2/v 



t=0 



We now replace the Vandermonde determinants in this expression with Vandermonde deter- 
minants of differential operators. Observe that 



A(w 2 ) = A[—) e EiU«?* 
dL 



Li=0 



and also 



A ( w ) . e £ti^ = A ( —) e ^ WiMi 
v 1 \dM J 



Mi=N 



This implies that we can compute the integral in the above formula as 



ps? / ■ • • / A{w)Aiw2) exp (£ ^ + 1) % 



2t \ dw\ dwk 



(rr/'' 



A I — ) A ( — 

dL/ V dM / ^ V.27ri j| w | =1 



1 f F ( ^ 2/^ \ ^t/j 



2t \ dwi (iwfc 



'J/.- w 2fc 



. , 2 . r 2t\ dw 
exp ( LjW + M jW + — 2 J ^ 



fcldet 



d i+j ~ 2 1 



fcx/c dL^dMi- 1 2m 



\w\=l 



exp ( Lw 2 + Mw + 



2t\ dw 



w 2 J w 2k 



L=0, 
M=N> 



where the last equality is by Lemma |6j Using the fact that 



d i+ i~ 2 1 



dL^dMi- 1 2m 



\w\=l 



exp ( Lw 2 + Mw + 



2t\ dw 



w 2 ) w 2k 



L=0, 
M=N 



1 / exp(Nw + 



it 



-dw 



2lX% J yjlk—2i j+3 
N 2k-2i- j+ 2 r Qxp ( w + 2tN^- 



2ixi 



w 



2fe-2i-i+3 



dw, 



we obtain a simplified formula for Mfc(USp(2iV), 2): 



kxk \ 2isi 



t=0 
2tN 2 - 
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2i-j+l 



4)*(*-W2fc !iV ^+5)/2^_ (Vdet f — / 6XP( 9 W + J 



•2u_ 



■dw 



t=o 



u=0 



where we have interchanged columns of the matrix to obtain the second line and set u = tN 2 
to obtain the third. The theorem follows. □ 



Proof of Theorem [H We can proceed in the same way as in the proof of Theorem [T], starting 
from Corollary [2j □ 
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Proof of Theorem^ This time we consider the third derivative. From fl6]) we have 



d 3k 



da\ ■ ■ ■ da\ \\ 



l<i,j<k\ w j 



w] - of) 



d 



2k 



OL\=— =Q!fc=0 



da\---da\ Ui<u< k ( w j ~ a f) 



a 1 =-=a k =0 



3 =l 3 \ j=l 3 J 

3=1 3 \ 3=1 J / 



t=0 



and the theorem follows from Corollary [3] in the same way. 



□ 



Proof of Theorem [^} We begin by proving a two- variable version of the recurrence relation. 
Let 



9m(x,y) :-- 



o XZ +72 



dz, 



2mJ lzl=l z^ 1 
Tk,e( x ^y) ■= det (g 2i - j+ e(x,y)) 

kxk 

for x, y G C, k > and £ G Z. Using Lemma UJ we will show that: 



d 2 ~ d ~ d ~ 

T k +w(x, y)Tk-i,i(x, y) = T k) e(x, y) T —^-T k/ {x, y) - T-T k/ (x, y)^-T k) e(x, y). (7) 

oxoy ox ay 



Let A k; i denote the matrix of T k /(x,y), i.e. 

93+1 



A 



9t 

92+e 



\ 92k-l+i 92k-2+i 



92-k+e \ 
94-k+e 

9k+t J 



Observe that 



d 

—g m (x,y) = <7 m _i(x,y,2t), 
d 



dy 



9m(x,y) = g m+2 (x,y,2t). 



We now compute the partial derivatives of Tk/(x, y). Expanding the derivative by columns, 
we obtain: 



TT" >k,t 
OX 



dg 1+e /dx 


gt 


• g2-k+e 






gi+e 


dgt/dx 


■ g2-k+i 


dg 


i+i/dx 


g2+e 


• c/4-k+e 


+ 




h+e 


dg 2+l /dx 


■ 9A-k+£ 


dg 2 k 


-i+i/dx 


fj2k-2+e ■ ■ 


(jk+e 




g2k~i+i 


dg 2k -i+e/dx ■ 


gk+i 






gi 


■ ■ dg 2 -k+t 


i/dx 








■ + 


g3+e 


g2+i 


■ ■ <9#4-fc+< 


i/dx 










£l2k-i+e 


g2k-2+e • 


dg k+e /dx 
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All terms but the last one in this sum vanish. Thus -§^7~ k ,e — det Ak+i ! e{ k ~^. 1 ) , where A^'"'^ 
denotes the matrix A with rows oi, . . . , a r and the columns b\, . . . , b s removed, as in Lemma 

m 

Similarly, expanding the partial derivative by rows, we have ^T k ,i = det A k +\£ ( fc+1 ) ■ 
Finally, jf^-Tke = det Ak+ieCl)- Thus Equation [7] is an immediate consequence of Lemma 

m y ' 

To return to our original functions g m (u) and %/(u), we observe that a simple change of 
variables gives g(x,y) = x m g m (x 2 y) for any x 7^ 0. Further, removing a factor of x 2l+e from 
each row i and x - - 7 from each column j in the determinant, we have 



If we set u = x 2 y, then by the chain rule we have = 2xy4-, S- = x 2 4~-, and so 

& 1 J nrr. < y rt.it 7 rhi ri.11. 7 



T k ,t(x,y) = x^ +kl T k/ (x 2 y). 

— = 2xv— 

dx * du 1 dy du 

%M d^y\ X TkAu) ) 

nk (k+l)+2kl+l 



X 



+k£ + 2) %MrUu) + 2uT k/ {u)T k "M 



and 

d ( Hh+i) hJ \ d 

Tx\ x 2 T vM)djH* * '« 

= a** 1 * 2 "* 1 ((^f^ + k?j %MTIM + 2u% t {u) 
Thus the theorem follows from (jS]) and (J7|). 

4. Numerical values 
Below are the first several values for the constant b k (USp(2N, 2)). See Theorem [TJ 

w 1 

b 2 



□ 



2-3 
19 



2 4 • 3 2 • 5 • 7 
487 

03 



2 7 ■ 3 5 ■ 5 2 ■ 7 - 11 
_ 59 • 197 

4 ~ 2 13 ■ 3 8 ■ 5 2 - 7 2 - 11 ■ 13 
174290791 

5 ~ 2 19 ■ 3 10 - 5 5 • 7 3 ■ ll 2 -13-17-19 
3373 • 1670407 

6 ~ 2 25 - 3 14 ■ 5 6 ■ 7 3 ■ ll 3 ■ 13 2 -17-19-23 

37-83-2203-3571457 

7 ~ 2 32 • 3 19 • 5 9 • 7 6 • ll 3 • 13 3 ■ 17 2 -19-23 

61 ■ 595351 ■ 11423948521 

2 42 . 3 23 . 5II . 7 7 . U 4 . 13 4 . 17 2 . 19 2 . 23 ■ 29 ■ 31 
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53 ■ 16646765854629827113 

" 2 53 • 3 29 • 5 13 • 7 9 • ll 5 • 13 4 • 17 3 • 19 3 • 23 2 -29-31 

47 ■ 1553 • 1787 ■ 73709 ■ 152825093 
2 62 . 334 . 517 . 710 . H 5 . 13 5 . 17 4 . 19 3 . 232 .29-31-37 



Below are several values for the constant bk{SO(2N),2). See Theorem [21 



61 = 1 
b 2 = ^- 

2-3-5 
2-13 

63 "3^7 
17-5987 

4 ~ 2 6 ■ 3 5 ■ 5 2 ■ 7 2 • 11 - 13 
_ 157-17519 

5 ~ 2 9 ■ 3 8 ■ 5 4 ■ 7 2 • 11 - 13 - 17 
22273664659 

6 ~ 2 15 ■ 3 12 • 5 5 • 7 3 • ll 2 • 13 2 -17-19 
_ 116228886131 

7 ~ 2 14 ■ 3 14 • 5 8 • 7 5 • ll 3 • 13 2 -17-19-23 
36774351481263481 

8 ~ 2 28 • 3 19 • 5 9 • 7 6 • ll 4 • 13 3 • 17 • 19 2 -23-29 
71 • 103 • 223 • 661 • 1069 • 134437 

9 - 2 34 . 325 . 511 . 76 . U 4 . 13 4 . 17 3 . 192 .23-29 
25171 • 7695491 • 57668937071891 

2 45 . 329 . 515 . 79 . H 5 . 13 5 . 17 3 . 19 3 . 23 2 -29-31-37 



We omit a table of values for the odd orthogonal case. Recall that by Theorems [Q and [31 
bk(0~(2N), 3) = 3 • 2 k ■ 6 fc (USp(2iV), 2) so that these values are given in terms of the above 
table for the symplectic case. 
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